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PREFACE 


Two  empirical  models  for  the  distribution  of  significant  wave  height  are 
given  in  Section  4.332  of  the  Shore  Protection  Manual.  A  similar  model  is 
presented  in  this  report.  The  model  is  based  on  a  three-parameter  Weibull 
distribution  function.  Parameters  in  the  model  are  evaluated  from  a  large 
sample  of  shallow-water  gage  data  at  Nags  Head,  North  Carolina.  The  model, 
which  more  closely  represents  available  data  than  either  of  the  previous 
models,  is  particularly  useful  for  statistical  prediction  of  extreme  signifi¬ 
cant  wave  heights  in  shallow  water  at  Nags  Head.  The  technique  is  applicable 
for  other  gage  sites.  This  work  was  carried  out  under  the  waves  and  coastal 
flooding  research  program  of  the  U.S.  Army  Coastal  Engineering  Research 
Center  (CERC) . 

This  report  was  prepared  by  Edward  F.  Thompson,  Hydraulic  Engineer,  under 
the  supervision  of  Dr.  C.L.  Vincent,  Chief,  Coastal  Oceanography  Branch.  The 
author  gratefully  acknowledges  Dr.  D.L.  Harris,  formerly  CERC  Senior  Scien¬ 
tist,  who  provided  valuable  comments  on  this  study,  and  J.  Peworchik  of  CERC 
who  processed  the  Nags  Head  data  for  the  study. 

Comments  on  this  publication  are  invited. 


Approved  for  publication  in  accordance  with  Public  Law  1966,  79th  Congress, 
approved  31  July  1945,  as  supplemented  by  Public  Law  1972,  88th  Congress, 
approved  7  November  1963. 


Colonel,  Corps  of  Engineers 
Commander  and  Director 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI)  UNITS  OF  MEASUREMENT 


f 


U.S.  customary  units 
metric  (SI)  units  as 

of  measurement 
follows : 

used 

in  this  report  can  be  converted  to 

Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0. 7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 

lO"3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angel) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins 

To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C  «  (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K  *■  (5/9)  (F  -32)  +  273.15. 


SYMBOLS  AND  DEFINITIONS 


d 

F[HSc  2  V] 


g 

H_ 


sc 


sc 


sc 


H 


sc  min 


T 

a 

a 


water  depth 

probability  that  the  ratio  Hg/Hg  =  Hgc  is  greater  than  or  equal 
to  a  specified  ratio  Hgc 

acceleration  due  to  gravity 

significant  wave  height 

mean  signifies  t  wave  height 

significant  wave  height  divided  by  mean  significant  height 
parameter  in  Weibull  distribution  function 
specified  value  of  Hgc 

minimum  expected  value  of  Hgc,  parameter  in  Weibull  distri¬ 
bution  function 

wave  period 

parameter  in  Weibull  distribution  function 
standard  deviation  of  significant  wave  height 
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A  MODEL  FOR  THE  DISTRIBUTION  FUNCTION 
FOR  SIGNIFICANT  WAVE  HEIGHT 

by 

Edward  F.  Thompson 

I .  INTRODUCTION 

The  long-term  distribution  of  significant  wave  height  at  a  site  can  be 
estimated  from  empirical  data  or  by  either  of  two  empirical  models  in  the 
Shore  Protection  Manual  (SPM) ,  Section  4.332  (U.S.  Army,  Corps  of  Engineers, 
Coastal  Engineering  Research  Center,  19 77) 1  .  Where  sufficient  data  are  avail¬ 
able,  direct  use  of  empirical  data  is  usually  preferable  to  either  SPM  model. 
However,  proper  estimation  of  the  long-term  distribution  of  significant  height 
requires  the  use  of  an  integral  number  (preferably  3  or  more)  of  reasonably 
complete  years  of  data.  This  requirement  is  often  difficult  to  meet  because 
of  intermittent  failure  to  obtain  observations  and  limited  number  of  years  of 
data  collection  at  a  site. 


It  is  often  convenient  to  model  an  observed  distribution  of  significant 
wave  height.  A  model  provides  a  simple  parameterization  of  the  observed  dis¬ 
tribution  as  well  as  a  systematic  method  for  extrapolating  to  probabilities 
beyond  the  data  (although  extrapolations  are  always  much  more  uncertain  than 
the  part  of  the  distribution  well  supported  by  data  because  of  long-term 
variability  in  storms  producing  extreme  wave  conditions).  Since  there  is  no 
compelling  physical  basis  for  favoring  any  particular  model,  models  are  chosen 
to  fit  observed  distributions  of  significant  height.  The  models  in  the  SPM 
were  proposed  as  a  tool  for  representing  the  distribution  of  the  highest  50  to 
80  percent  of  observed  significant  heights.  The  model  presented  in  equation 
4-6  of  the  SPM  is  a  two-parameter  modified  exponential  distribution  which  is 
further  simplified  in  equation  4-9  of  the  SPM  to  a  one-parameter  distribution. 


The  model  presented  in  this  report  is  based  on  a  three-parameter  Weibull 
distribution  function.  The  three-parameter  model  can  better  fit  observations 
than  either  the  one-  or  two-parameter  models  in  the  SPM.  Parameters  are  eval¬ 
uated  to  optimize  the  fit  to  empirical  data  from  a  gage  at  Nags  Head,  North 
Carolina.  The  model  is  formulated  in  dimensionless  terms  so  that  the  effect 
of  mean  significant  wave  height  level  is  removed.  The  advantages  of  using 
dimensionless  terms  are  that  more  complete  use  is  made  of  available  data  and 
that  general  characteristics  of  the  distribution  of  significant  height  in 
addition  to  the  mean  can  be  readily  examined.  The  dimensionless  distribution 
function  may  be  relatively  invariant  compared  to  mean  significant  height  vari¬ 
ations  along  a  short  section  of  coast.  Hence,  the  model  presented  in  this 
report  is  believed  to  provide  a  more  general  representation  than  the  models 
in  the  SPM. 


!U.S.  ARMY,  CORPS  OF  ENGINEERS,  COASTAL  ENGINEERING  RESEARCH  CENTER,  Shore 
Protection  Manual,  3d  ed.,  Vols.  I,  II,  and  III,  Stock  No.  008-022-00113-1, 
U.S.  Government  Printing  Office,  Washington,  D.C.,  1977,  1,262  pp. 


7 


II.  WEIBULL  DISTRIBUTION  FUNCTION 


A  general  model  which  can  be  used  to  approximate  the  empirical  distribution 
of  significant  wave  height  is 


FfH  > 
L  sc 


(1) 


where 


H 

H 

F 

H 

H 


sc 

sc 

[H-c* 

sc  min 


significant  wave  height  divided  by  mean  significant  height 
specified  value  of  Hgc 

probability  that  Hgc  is  greater  than  or  equal  to  a  speci¬ 
fied  ratio  Hg<=> 

minimum  expected  value  of  Hgc 

other  parameters  in  distribution  function. 


Equation  (1)  is  a  form  of  the  Weibull  distribution  function  with  three  param- 
eters  (IIgc  min,  Hsc,  and  a). 


III.  COMPILATION  OF  EMPIRICAL  DISTRIBUTION 

The  parameters  in  equation  (1)  must  be  evaluated  for  each  site  by  opti¬ 
mizing  the  agreement  between  equation  (1)  and  the  empirical,  dimensionless 
distribution  of  significant  height  at  the  site  using  the  following  procedure: 

(a)  Assemble  all  significant  heights  obtained  by  reliable,  consist¬ 
ent  analysis  methods  at  a  particular  site; 

(b)  delete  significant  heights  from  months  in  which  more  than  50 
percent  of  the  possible  observations  are  missing; 

(c)  compute  mean  significant  height  for  each  remaining  month; 

(d)  divide  each  significant  height  by  the  appropriate  monthly  mean 
significant  height;  and 

(e)  combine  the  dimensionless  significant  heights  in  step  (d)  from 
all  months  into  one  distribution. 

The  implicit  assumption  in  step  (e)  is  that  monthly  variations  in  wave  condi¬ 
tions  can  be  completely  represented  by  variations  in  monthly  mean  significant 
height  but  that  variations  relative  to  the  mean  are  consistent  from  month  to 
month.  The  assumption  may  not  be  valid  at  sites  strongly  affected  by  hurri¬ 
cane  waves  unless  hurricane  waves  are  treated  separately. 
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IV.  APPLICATION  TO  SHALLOW-WATER  GAGE  SITE 
AT  NAGS  HEAD,  NORTH  CAROLINA 


Wave  data  used  to  test  the  model  were  obtained  from  a  pier-mounted  staff 
gage  in  a  16-foot  water  depth  at  Nags  Head,  North  Carolina  (see  Thompson, 
1977)2.  Digital  records  were  collected  and  analyzed  at  6-hour  intervals,  with 
numerous  interruptions,  from  December  1968  to  March  1978.  Significant  wave 
heights  from  each  of  54  relatively  complete  months  were  processed,  as  discussed 
previously,  to  form  one  dimensionless  distribution  containing  5,220  observa¬ 
tions  (see  Fig.).  The  empirical  distribution  extends  to  an  exceedance  percent¬ 
age  of  0.02.  Many  cases  at  low  exceedance  percentages  may  be  affected  by  the 
limited  water  depth  at  the  Nags  Head  site,  as  indicated  by  the  hatched  area  in 
the  Figure. 
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Figure.  Distribution  of  significant  wave  height. 
Nags  Head,  North  Carolina. 


2TH0MPS0N,  E.F. ,  "Wave  Climate  at  Selected  Locations  Along  U.S.  Coasts," 
TR  77-1,  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center, 
Fort  Belvoir,  Va.,  Jan.  1977. 
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Values  for  the  parameters  in  equation  (1)  were  estimated  from  the  empirical 
distribution  of  significant  heights  at  Nags  Head  by  the  method  described  in  the 
Appendix  to  give 


r  -  -/hsc  -  0.198V'65 

LHsc  -  Hsc]  =  e  \  0.885  ' 

(2) 

Equation  (2)  is  shown  in  the  Figure.  This  distribution  function  fits  the 
empirical  distribution  at  Nags  Head  better  than  the  comparable  models  in  the 
SPM.  Equation  (2)  can  be  rearranged  to  give 

jo. 606  Jt.n (-  Hn  fIh^  >  Hc  1)-  0.122 
Hsc  =  ej  V  LSC  SCJ/  1+0.198  (3) 

The  assumption  that  variations  relative  to  monthly  mean  significant  height 
are  consistent  from  month  to  month  was  tested  at  Nags  Head  by  comparing  empir¬ 
ical  distributions  of  significant  height  by  season.  Distributions  for  fall 
(September  to  November) ,  winter  (December  to  February) ,  and  spring  (March  to 
May)  are  comparable  to  the  empirical  distribution  in  the  Figure.  However,  the 
distribution  for  summer  (June  to  August)  indicates  higher  values  of  Hgc  than 
the  empirical  distribution  in  the  Figure  at  probabilities  below  about  10  per¬ 
cent.  This  discrepancy  is  due  to  exceptionally  low  monthly  mean  significant 
heights  in  summer  and,  in  many  cases,  to  nearshore  hurricanes.  Equation  (2) 
seems  to  be  satisfactory  for  estimating  the  annual  or  nonsummer  distribution 
of  significant  height  at  Nags  Head,  but  the  summer  distribution  must  include 
special  consideration  of  hurricane-generated  waves  at  this  site. 


V.  EXAMPLE  PROBLEMS 

***************  EXAMPLE  PROBLEM  1*********  ****** 

GIVEN :  Mean  annual  significant  height  of  approximately  Hs  =  3.0  feet  (0.91 
meter)  at  Nags  Head,  North  Carolina  (see  Table). 

FIND: 


(a)  The  significant  height  which  is  equaled  or  exceeded  during  6 
hours  every  year. 

(b)  The  significant  height  which  is  equaled  or  exceeded  during  1 
hour  every  year. 

SOLUTION: 

(a)  The  exceedance  percentage  (F[HSC  >  Hsc]  x  100  percent)  is 
1 


365  x  24 
6 


x  100  =  0.0685  percent 


Table.  Monthly  and  annual  significant  wave  height  statist! c«s ,  Nags  Head, 
North  Carolina 


No .  of 

No .  of 

Month 

Year 

Ohs ns . 

«s 

a 

Month 

Year 

Obsns . 

Hs 

0 

(ft) 

(ft) 

(ft) 

(ft) 

12 

1968 

96 

2.62 

1.34 

04 

1974 

73 

2.50 

1.19 

01 

1969 

112 

3.04 

1.59 

09 

1974 

75 

2.55 

0.98 

02 

1969 

103 

4.12 

2.01 

12 

1974 

118 

3.09 

1.59 

03 

1969 

106 

3.66 

1.85 

01 

1975 

93 

3.  38 

1.58 

04 

1969 

90 

2.43 

1.27 

03 

1975 

63 

2.97 

1.72 

05 

1969 

85 

2.32 

1.08 

08 

1975 

63 

1.86 

0.86 

07 

1969 

1 12 

2.01 

1.02 

09 

1975 

86 

2.99 

1.10 

08 

1969 

94 

2.29 

1.07 

10 

1975 

92 

3. 35 

1.40 

09 

1969 

105 

3.36 

1.74 

1 1 

1975 

90 

3.02 

1.49 

Annual 

Dec.  1968-Oct.  1969 

1,006 

3.00 

1.72 

02 

1976 

87 

2.30 

1.09 

09 

1971 

117 

3.44 

1.78 

03 

1976 

99 

2.71 

1.13 

10 

1971 

117 

3.34 

2.02 

04 

1976 

72 

2.58 

1.33 

11 

1971 

78 

3.88 

1.71 

06 

1976 

65 

2.06 

1.01 

12 

1971 

120 

3.39 

1.94 

10 

1976 

113 

2.85 

1.14 

01 

1972 

82 

2.38 

0.90 

11 

1976 

113 

2.21 

1.19 

02 

1972 

116 

3.81 

2.03 

12 

1974 

109 

2.63 

1.26 

03 

1972 

123 

2.93 

1.21 

01 

1977 

89 

2.90 

1.18 

04 

1972 

110 

3. 10 

1.67 

02 

1977 

98 

2.21 

1.18 

05 

1972 

120 

3. 19 

1.53 

03 

1977 

118 

2.34 

0.93 

06 

1972 

101 

2.15 

0.98 

04 

1977 

93 

1.97 

0.78 

08 

1972 

88 

2.24 

1.33 

05 

1977 

82 

1.60 

0.72 

Annual 

Sept.  1971-Aug.  1972 

1,173 

3.  10 

1.71 

06 

1977 

91 

1.68 

0.82 

09 

1972 

106 

3.06 

2.03 

08 

1977 

88 

1.40 

0.59 

10 

1972 

109 

3.77 

1.84 

12 

1977 

93 

3.20 

1.40 

11 

1972 

96 

4.58 

1.40 

01 

1978 

100 

3.26 

1.21 

12 

1972 

97 

2.96 

1.61 

03 

1978 

82 

3.48 

1.23 

01 

1973 

97 

3.36 

1.61 

02 

1973 

92 

4.43 

2.09 

03 

1973 

114 

4.33 

2.28 

05 

1973 

89 

2.38 

1.00 

Annual 

Sept.  1972-May  1973 

854 

3.55 

1.93 

1  1 


From  the  Figure  or  equation  (3) 


Hsc  =  =  3.13, 


Hg  =  9.4  feet  (2.9  meters). 


(b)  The  exceedance  percentage  is 

365  x  24  X  =  ^.01 14  percent. 
From  the  Figure  or  equation  (3), 


Hsc  =  f  =  3‘57’ 
Hs 


Hs  =  10.7  feet  (3.3  meters). 


Check  to  see  if  H_  exceeds  depth-limited  height. 


d_ 

H 


16 

10.7 


=  1.5 


From  Figure  2-66  in  the  SPM,  depth-limited  breaking  may  be  possible  if 


gT2 


>  0.0172 


This  condition  corresponds  to 


T 


»s 

0.0172  g 


10.7 


0.0172  x  32.2 


4.4  seconds. 


Since  a  period  of  4.4  seconds  or  less  is  unreasonably  short  for  a  10.7-foot- 
high  wave  at  this  site  (see  Thompson,  19 77) 3 ,  depth-limited  breaking  is  not 
expected  to  be  a  consideration  in  this  example. 

***************  EXAMPLE  PROBLEM  2 *************** 

GIVEN :  Mean  significant  height  of  approximately  Hs  =  3.4  feet  (1.0  meter)  in 
February  at  Nags  Head,  North  Carolina  (see  Table). 

FIND:  The  significant  height  which  is  equaled  or  exceeded  during  6  hours  every 
February . 

SOLUTION :  The  exceedance  percentage  is 


1 

— -r - r  x 

28  x  4 


100  =  0.89  percent 


3TH0MPS0N,  E.F.,  op.  cit.,  p.  9. 
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From  the  Figure  or  equation  (3) 


*sc 


Hs  =  8.4  feet  (2.6  meters). 

VI.  SUMMARY 

A  three-parameter  model  for  the  distribution  of  significant  wave  height  is 
given.  A  procedure  for  using  available  data  from  a  site  to  compile  a  dimension¬ 
less  distribution  of  significant  height  and  to  estimate  parameters  in  the  model 
is  presented.  The  procedure  extends  the  use  of  available  data  and  leads  to  a 
model  which  more  closely  follows  the  data  than  procedures  given  in  the  SPM.  The 
procedure  is  applied  to  shallow-water  gage  data  from  Nags  Head,  North  Carolina. 
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APPENDIX 


METHOD  FOR  ESTIMATING  PARAMETERS  IN  THE  WEI BULL  DISTRIBUTION  FUNCTION 


The  Weibull  distribution  function,  equation  (1),  can  be  transformed  into  a 
form  suitable  for  linear  regression  analysis.  First,  the  natural  logarithm  of 
equation  (1)  is  taken 


£n  F 


(A- 1) 


Both  sides  of  equation  (A— 1 )  are  multiplied  by  -1,  and  again  natural  logarithms 
are  taken 


£n  (-  £n  F)  =  in 


H  -  H  \a 

sc  sc  min) 

*  Hsc  / 


Equation  (A-2)  can  be  rewritten  as 

(bsc  “  ^sc  min 


(A-2) 


)  =  Hsc  +  a 


in  (-  in  F) 


(A- 3) 


Equation  (A-3)  is  in  the  form 


Y  =  a  +  b  X, 


(A-4) 


where 

Y  -  5,ri(nsc  -  Hgc 
a  =  in  Hsc 

b  -  i 

X  =  in  (-  in  F)  . 


An  initial  value  of  the  parameter  Hsc  was  obtained  from  Table  1  of 

Thompson  and  Harris  (1972) 4  as  the  "minimum  significant  height"  divided  by  the 
observed  mean  significant  height.  The  value  for  Nags  Head  was  0.31.  Alterna¬ 
tively,  Hsc  m^n  could  be  estimated  initially  as  0.38  from  equation  4-8  in  the 
SPM.  The  estimated  value  of  Hsc  m^n  and  empirical  tabulation  of  F  as  a 
function  of  Hsc  are  used  to  compute  a  table  of  X  and  Y  values.  Linear 
regression  analysis  Is  then  used  to  estimate  optimum  values  of  a  and  b  in 


^THOMPSON,  E.F.,  and  HARRIS,  D.L.,  "A  Wave  Climatology  for  U.S.  Coastal 
Waters,"  Proceedings  of  Offshore  Technology  Conference ,  May  1972,  pp.  675-688 
(also  Reprint  1-72,  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research 
Center,  Fort  Belvoir,  Va. ,  NTIS  AD  746  365). 
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equation  (A-4).  Values  for  the  parameters  Hsc  and  a  in  equation  (1)  are 
easily  calculated  from  a  and  b  by  using  the  relationships  defined  for 
equation  (A-4). 

The  distribution  function  given  by  equation  (1)  with  the  estimated  para¬ 
meters  Hsc  min,  Hsc,  and  a  is  compared  with  the  empirical  distribution. 

A  new  value  of  Hsc  min  is  estimated  to  attempt  a  better  fit  to  the  empirical 
distribution.  The  new  Hsc  m^n  is  used  to  compute  a  new  table  of  X  and  Y 
values  which  is  then  used  to  estimate  new  values  of  Hgc  and  a  as  before. 
The  process  is  continued  until  an  Hsc  m^n  has  been  found  which  leads  to  a 
satisfactory  model.  The  parameters  for  Nags  Head  in  equation  (2)  are  con¬ 
sidered  satisfactory  because  they  specify  a  distribution  function  which  fits 
the  lower  99  percent  of  the  empirical  distribution  reasonably  well  but  is  con¬ 
servatively  high  in  comparison  to  the  highest  1  percent  of  the  empirical  dis¬ 
tribution.  Since  the  highest  part  of  the  empirical  distribution  is  based  on 
the  smallest  number  of  observations,  it  is  the  least  well-established  part  of 
the  empirical  distribution,  and  conservatism  is  desirable  in  a  model  for  engi¬ 
neering  use. 


i 


